BCS-BEC crossover in a random external potential 
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We investigate the ground state properties of a disordered superfluid Fermi gas across the BCS- 
BEC (Bose Einstein condensate) crossover. We show that, for weak disorder, both the depletion of 
the condensate fraction of pairs and the normal fluid density exhibit a nonmonotonic behavior as a 
function of the interaction parameter 1/fcpa, reaching their minimum value near unitarity. We find 
that, moving away from the weak coupling BCS regime, Anderson's theorem ceases to apply and 
the superfluid order parameter is more and more affected by the random potential. 

PACS numbers: 



Ultracold Fermi gases near a Feshbach resonance have 
been the subject of intense experimental and theoretical 
investigations over the last five years [if. By tuning the 
inverse s-wave scattering length 1/a continuously from 
-oo to +oo, the ground state of these systems trans- 
forms from a weak-coupling BCS superfluid into a Bose- 
Einstein condensate of diatomic molecules. 

An interesting and general problem is to explore the 
BCS-BEC crossover in the presence of a random external 
potential. Optical disorder can be introduced in ultracold 
gases in a highly controllable way. For instance, in recent 
experiments with elongated atomic BEC 0, d, 0] , speckle 
patterns, created by a laser beam passing through a dif- 
fusive plate, were shown to suppress the ID expansion 
of the gas. A characteristic feature of these random po- 
tentials is that they are slowly varying in space, with a 
typical length scale of the order of a few /im. In order to 
produce a real white noise disorder, which is ubiquitous in 
condensed matter theory, it was recently suggested [5| to 
(randomly) trap atoms of a different species at the nodes 
of a tight 3D optical lattice. If the optical polarizability 
of the fermions is small, the latter will not experience 
the periodic potential but only the interaction with the 
impurity atoms. 

To date, the effects of impurities on a superfluid Fermi 
gas near unitarity (1/a = 0) are unknown, even for weak 
disorder. For s-wave BCS superconductors, Anderson's 
theorem @ states that the (disorder-averaged) order pa- 
rameter is unaffected by nonmagnetic impurities, as long 
as localization effects can be neglected. The superfluid 
density, measuring the dynamical response of the system 
to a superfluid flow, is instead reduced by impurities 0. 
For Bose superfluids, the depletions of the condensate 
fraction and of the superfluid density induced by a weak 
disorder have been calculated in Refs [1, [§| . 

In this Letter we provide the first theoretical study of 
a disordered superfluid Fermi gas across the BCS-BEC 
crossover, at zero temperature. We show that the de- 
pletion of the condensate fraction of pairs and the nor- 
mal fluid density, induced by a weak disorder, exhibit 
nonmonotonic behavior as a function of the interaction 
strength, with a pronounced minimum near unitarity. 



These results are therefore consistent with the general 
expectation that superfluidity is more robust against sup- 
pressing mechanisms in the crossover regime. Nonmono- 
tonic behaviors across the resonance have recently been 
observed [Io| in critical velocities for superfluid flow, in 
agreement with earlier theoretical predictions 11 , 13] • 

Our work is based on the Nozieres and Schmitt-Rink 
(NSR) theory [13] of superconducting fluctuations, ex- 
tended to the broken symmetry state 14{. In particu- 
lar, we will make use of techniques developed in Ref . fTHI] . 
where the NSR theory has been applied to calculate the 
condensate fraction and the normal fluid density of clean 
Fermi gases across the BCS-BEC crossover. 

In the following we assume that the random potential 
V(r) originates from the scattering of fermions against 
impurity atoms and is given by V(r) = J2j 9d5{r — R,-), 
where gd is the fermion-impurity coupling constant and 
Rj are the static positions of the impurities. The cor- 
responding correlation function takes the white noise 
form (V(—q)V(q)) = ftS^fiK, where q = (q, iuj m ) and 
k = nig\, rii being the concentration of impurities. 

The two-component Fermi gas is described by the 
imaginary time action S = Jq dr J" dr{-0 cr c?T-'(/; cr + H}, 
where (3 in the inverse temperature and 



H = {K + V) V> 



(1) 



Here K = — V 2 /2m — /i, /.i is the chemical potential, g is 
the fermion-fermion interaction and we set h = Volume = 
k B = 1. 

Following Ref. [IBj , we apply the Hubbard-Stratonovich 
transformation to the Hamiltonian ([T]), by introducing 
a new bosonic field A(x) which couples to tp^l- I n_ 
tegrating out the fermionic fields, the partition func- 
tion of the system takes the path integral form Z = 
J D[A, A]e~ s eft , where 



S, 







dr / dr 



\A[x)l 
9 



Tr\n[-f3G~ 1 {x)} 



. (2) 

is the effective action written in terms of the inverse single 
particle Green's function 

G~\x) = -d T - (K + V(r))r 3 + A(x)t+ + A(x)t- . (3) 



2 



Here x — (r, r), Tj are the Pauli matrices and t = 
(n±ir 2 )/2. 

To procede further, we assume that the main contri- 
bution to the partition function comes from small fluc- 
tuations SA(x) = A(x) — A around the uniform BCS 
pairing field A. Hence, we expand the effective action 
([2]) up to quadratic order in the bosonic fields SA and 
in the random potential V. To this purpose, the Green's 
function |3]) is conveniently written as G _1 = Gq 1 + S, 
where Gq 1 = — d T — Kt% + Ati is the inverse BCS 
Green's function in the absence of disorder, and S = 
-V(t)t 3 + SA(x)t + + 5A(x)t~ . 

Averaging over the positions of the impurity atoms, 
we obtain S cS = S F + A(S A(0) + M(0)) + S B , where 
S F = {3A 2 /g - /T 1 E fe TrH-GoW + V(0)r 3 ]_+ pQ F 
corresponds to the mean field action. Here G 1 (fc) — 
iuj n — ^k7"3 + Ati is written in momentum representation, 
with £k = |k| 2 /2m — [i and 



kq 



(4) 



where G = G (fc) and Gq = G (k + q). Further- 
more, A = A(l/g- X; k tanh(/3£: k /2)/2i; k ), where £ k = 
+ A 2 ) 1//2 is the energy spectrum of single particle 
excitations, and 

Sb = \H [r) f Mv + V(-q)W^r) + V(q)Wrf] (5) 



corresponds to the gaussian action for the bosonic fluc- 
tuations. Here rf = [8A(q], SA(— q)] and M is a 2x2 
symmetric matrix whose elements are given by 

M n (q) = - + J2 G °'™ G 'w M M =E G o. 12G M2' 

" fc k 

(6) 

and M22(q) = Mn(—q) (explicit expressions can be 
found in Ref.[H|). In Eq.©, the doublet 



W(q) 



En fit ft 
k «J0,12<J0,11 _ ^0,22(^0,12 

En fit n t~*i 

u Ltq, 11^0,12 — ^0, 12^0,22 



(7) 



couples disorder to the bosonic fluctuations and induces 
processes in which pairs are scattered in or out of the con- 
densate. At zero temperature, W\ = W2 — J^k A(£k + 
£k +q )/2£ k £ k+q (£ k + £ k+q ), for q = (q, 0). 

By assumption, the mean field action Sp is an ex- 
tremum of S e g, hence the coefficient A must be set 
equal to zero (A = 0). This yields the usual BCS gap 
equation 1/g — tanh(/3_B k /2)/2i? k , as in the ab- 
sence of impurities. This equation is ultraviolet diver- 
gent and is cured [14] by replacing the bare coupling 
constant g with the scattering length a according to 
1/g = —m/Aira + l/2eit) with e k = k 2 /2w. At zero 



temperature, one finds 



m 
Aira 



E 



1 

2E^ 



1 

2^ 



which should be seen as a self-consistent relation between 
A and [i. 

In the NSR theory, the second constraint is provided by 
the number equation n = —dtl/d/i, where Q = — ]xiZ//3 
is the thermodynamic potential. To calculate it, we sub- 
stitute S e ff = Sf + Sb in the partition function Z and 
integrate over the Bosc fields. At zero temperature, this 
gives f2 = fi_F + where £If = Sf/ (3 is the mean field 
fermionic term and 



SIr = lim 

/3^oo 2/3 



-^lndetM-J ^ W^M^W (9) 



q,w m — 



gives the contribution from bosonic fluctuations. 

The integration over q in Eq.f4|) is also ultraviolet di- 
vergent, due to the zero range nature of the fermion- 
impurity potential. This problem is solved by a renor- 
malization of the related coupling constant g& — > g& + 
9d !/ e q- Equation (gj) then yields tt d F = kA 2 to 3 /47t 2 , 
which is independent of /1. As a result, the number equa- 
tion takes the form 



k 



1 - 



(10) 



Differently from Eq.©, Eq.fTU]) depends explicitly on 
disorder through the contribution © from bosonic fluc- 
tuations. From Eqs (|SJ) and (JTUJ) we obtain A and /jasa 
function of the atom density n — kp/3ir 2 , the scattering 
length a and the disorder strength k (up to linear terms). 

In the BCS limit, corresponding to 1/kpa — * —00, 
the bosonic contribution in Eg. (TIT)]) is negligible and 
we find fj, — ep — hp /2m while Eq.© gives A = 
8e~ 2 ep exp(— 7r/2fcp|a|), as in the absence of disorder. 
This means that for weak attraction the superfluid order 
parameter is unaffected by the disorder, in agreement 
with Anderson's theorem [6|. 

In the opposite strong coupling regime, correspond- 
ing to l/kpa — * +00, the gap equation ([8]) yields /1 — 
— l/2ma 2 , which is half of the binding energy. To find 
A, we first differentiate Eq.© with respect to /1 and 
then expand the obtained result in powers of q, as- 
suming |q| 2 /2m, iuo m <C Taking into account that 
Wi(q = 0) = m 2 aA/47r, from Eq.([10]) we obtain 



n 
2 



m 2 A 2 a 



E 

q 



fJ-M 



2uj c 



-4k 



m 2 A 2 a/87r 



q 



2^ A f) 2 ' 
(11) 



where e 



M 



2 /4m, 



-Ml 



- q - q /4m, W q - 
trum of the bosonic excitations 



2[/,m) is the spec- 

In Eq.(|n]), Hm = 

A 2 a 2 m/2 — 2-Kan/m is the effective chemical potential 
for the composite bosons corresponding to a molecular 
scattering length cim — 2a [l|| . 

Equation has a direct physical meaning: the term 
(g) (m 2 A 2 a/87r) corresponds to the fraction of condensed 
molecules whereas the second term yields the Bogoliubov 
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FIG. 1: Quantum depletion a — ao = — (Km 2 /fcf)Q(l/fei?a) 
of the condensate fraction of pairs a induced by disorder: the 
function Q(l/fc_Fa) is shown with the solid line. The asymp- 
totic behavior at strong coupling is shown with the dashed 
line. Notice that Q saturates to 3/2n 2 at weak coupling (ar- 
row). 



quantum depletion due to the molecule-molecule inter- 
action, as shown in Ref. [HI]. The last term in Eq. (fTT|) 
gives the density of molecules pushed out of the conden- 
sate by the random potential. In particular, the fac- 
tor 4 accounts for the fact that the effective disorder 
potential Vm seen by the molecule is Vm(t) — 2F(r), 
hence the corresponding correlation function is given by 

(Mq)Viif(-q)}=4/s. 

From Eq. lfTTj) we also find A — Ao = —2nm/na : where 
Ao is the order parameter in the absence of disorder, 
showing that Anderson's theorem does not apply in the 
BEC limit. In particular, the correction A — Ao diverges 
for a — > + . This signals that in the absence of repulsive 
interactions, the uniform Bose gas of molecules is unsta- 
ble against collapse into a localized orbital of the disorder 
potential @- 

Condensate fraction of pairs. This is defined as the 
dimensionless ratio a — 2n c /n, where n c = J \G%2(r, t — 
0)| 2 dr is the density of condensed pairs corresponding to 
the normalization of the Cooper pair wavefunction. In 
the absence of disorder, the condensate fraction is an in- 
creasing function of the interaction strength It 
has been calculated within the NSR approach in Ref. 15 1 
and compared with previous mean field [13] and Quan- 
tum Monte Carlo [l8( predictions. 

We notice that disorder enters the formula G" 1 = 
Gq" 1 + X in two different ways: i) directly, through the 
"self energy" S, and ii) in the mean field propagator Go, 
through the shifts 5fi = p — po, SA = A — A of the vari- 
ational parameters, as calculated from Eqs dHJ and (flDj). 
Hence, the change in the density of condensed pairs is 
given by 



where n™ f = ^ k A 2 /4E 2 . is the mean field prediction for 
the condensate fraction and 8G d (k) corresponds to the 
disorder correction of the Green's function obtained by 
expanding the latter in powers of S and retaining up to 
second order terms. 

We have calculated the rhs of Eq. (fT2")) numerically. We 
find a — a = — (urn 2 /kp)Q(\/kpa), where the crossover 
function Q(l/kpa) is plotted in FigfTJwith the solid line. 
We see that the depletion of the condensate fraction ex- 
hibits a pronounced minimum at l/kpa « —0.9. 

In the deep BCS limit the leading contribution to the 
depletion (TT2|) comes from 5Gf 2 - Neglecting bosonic 
fluctuations, we have 5G d (k, iuj) = Got 3 GoT 3 Go. 
From Eq. (fl"2")) we then find n c — n c0 = —nm 2 k F /ATr 4 
corresponding to Q — 3/2ir 2 , as shown by the arrow 
in FigfTJ In the clean limit the condensate fraction 
n C Q = mkpA/S-n vanishes exponentially as l/kpa — > 
—00, so the effects of disorder are clearly important in 
this regime, despite the fact that the order parameter it- 
self is not affected by impurities. This interesting result 
follows by noticing that, in the presence of disorder, the 
zeroth-order Green's function Go,i2(r,T = 0) in coordi- 
nate space is damped by the factor exp(— r/2£), where 
£ — 7r / Km 2 is the mean free path between collisions 0]. 
Hence the order parameter A = gGi2{r = 0,0) is the 
same as in pure superfluids, while the condensate fraction 
n c = J |G„,i2(r,0)| a e-^dr ~ / \G , 12 (t, 0)| 2 (1 - r/£)dv 
is severely suppressed by disorder. Here the condition of 
weak disorder (n c o — n c )/n C Q <C 1 yields Krnfcr/A <C 1. 

In the BEC regime, the leading contribution to 
the depletion (|12p comes from the disorder induced 
change 5 A of the order parameter and we recover n c — 
n c o = — A«Am 3 /27r 2 , in agreement with Eq. pTjl and 
Ref. [8(. The corresponding asymptotic behavior Q = 
(12/n) 1 / 2 / \fk~F~a is shown in FigQ]with the dashed line. 
We see that the molecular gas is fully depleted if the 
scattering length vanishes. Notice however that weak 
disorder implies Km /4 7 v/ 2 « 1, so one cannot take 
the limit a — > + for fixed disorder k. 

Normal fluid density. This is defined as p n = n — p s , 
where p s is the superfluid density. The latter is related 
to the lowest order change f2[Q 2 ] — f2[0] = p s mv 2 /2 in 
the thermodynamic potential of the gas in the presence 
of a supercurrent with velocity v s = Q z /2m along, say, 
the z-axis. At zero temperature and in the absence of 
disorder, the superfluid density coincides with the total 
density. Therefore p„ = —\imQ^QAmd 2 Vl d [Q z ]/dQ 2 , 
where f2 d [(5 z ] is the contribution to the thermodynamic 
potential coming from disorder. 

To apply this formula, we notice that the current- 
carrying state is related to the equilibrium state by the 
gauge transformation A(r) — > A(r)e l ® zZ . Accordingly, 
the BCS Green's function Go maps into Go, where [l5| 
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FIG. 2: Zero temperature normal density p n /n = 
{nm 2 /k,F)F{l/kFa) induced by disorder: the function 
F(l/kFa) is plotted with the solid line. The asymptotic be- 
haviors at weak and strong coupling are shown with dashed 
lines. 



Repeating the calculations from Eq.Q to Eq.@ with 
Go replaced by Go, we obtain 

n d lQ*} = ~£>[Go7- 3 G&t 3 ] -^W^M-iff, (14) 

q,k q 

where Go = Go(fc) and G = Go(k + q), with q = (q, 0). 
In Eq. (fT4"|) , M and W are given by the expressions © and 
(O, respectively, after replacing Go, G with Go, G' . 

Differentiating twice Eq.([T3|l with respect to Q z and 
performing the numerical integrations, we find p n /n = 
(Km 2 /kF )F(l/kpa), where the function F (l/kpa) is 
plotted in Figf2] with the solid line. We see that also 
the normal fluid density exhibits a nonmonotonic be- 
havior, reaching its minimum value near unitarity, at 
l/kpa w -0.5. 

In the BCS regime, the leading contribution to the 
potential (fT4")) comes from the first term in the rhs 
and we find p n — K,mk F /24:ir 2 A, in agreement with 
Ref.Q. The corresponding asymptotic behavior F — 
(e 2 /32) exp(7r/2fcf|a|) is shown in Fig[2]with the dashed 
line. In this limit the healing length of the gas is defined 
as £bcs = kp/mA and we obtain p n /n — it£bcs /8^, 
showing that the normal density is proportional to the 
healing length for £bcs <C L 

In the opposite BEC limit (l/kpa ^> 1), the major 
contribution in Eq. (|14|) comes from the bosonic fluctua- 
tions. By expanding W , M in powers of q as done be- 
fore Eq. (fTTj) . we find p n = /t4m 3 A/37r 2 , which is con- 
sistent with Ref. [8j. The corresponding asymptotic be- 
havior F = (8/v37r) / ^kpa is plotted in Figf2]with the 
dashed line. 

In conclusion, we have discussed the ground state prop- 
erties of a disordered superfluid Fermi gas across the 



BCS-BEC crossover. We have shown that superfluidity 
and condensation of pairs are less affected by disorder 
near unitarity. Remarkably, in the BCS limit an expo- 
nentially weak disorder strongly suppresses the conden- 
sate fraction while the order parameter itself is insensitive 
to impurities, in agreement with Anderson's theorem. 

An interesting direction for future work is to generalize 
the above theory at finite temperatures and to study the 
effect of impurities on the superfluid transition tempera- 
ture T c . For BCS superfluids, the thermodynamics of the 
gas (and therefore T c ) is insensitive to impurities, due to 
Anderson's theorem. In contrast, for weakly interacting 
Bose condensates, the transition temperature is affected 
by disorder [l9j ]. 
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